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Abstract 

We investigate a phase transition of the 0(N) invariant scalar model using the 
auxiliary mass method. We determine the critical exponent (3 by calculating an effective 
potential below the critical temperature. This work follows that of a previous paper. % 
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Phase transitions at finite temperature are important phenomena in particle physics, 
cosmology and condensed matter physics. For example, the QGP phase is produced in 
heavy ion collisions. & Some phase transitions occured in the early universe. & The electro- 
weak phase transition in particular plays an important role in the electro-weak baryogenesis 
scenario^ and gives some constraints to models of elementary particle physics, i^'&'i-' We 
also see a great number of phase transitions in condensed matter physics. In the present 
paper, we investigate an 0(N) invariant scalar model which corresponds to many condensed 
matter systems, for example alloys, superfluids, and binary liquids. & 

To investigate such phase transitions, we can use finite temperature field theory, which 
is based only on a statistical principle. However, we often have an infrared divergence and 
cannot obtain reliable results using perturbation theory at finite temperature. & To over- 
come this problem, we used the auxiliary-mass method, and calculated an effective 
potential and critical exponents of the 0(N) invariant scalar model above the critical tem- 
perature Tc in a previous paper. We did not investigate at a temperature below Tc for 
two reasons, numerical instability and the lack of computer power. In this work we have 
overcome these problems, and we calculate an effective potential and critical exponents of 
the 0(N) invariant scalar model below the critical temperature. 

We explain the idea of the auxililary-mass method. Since we can calculate a reliable 
effective potential for temperatures T <C y using perturbation theory, B first we assume the 
mass as m ~ T and calculate an effective potential. This potential is reliable if the coupling 
constant. A, is small. We next extrapolate the effective potential to that of a true mass 
using a non-perturbative evolution equation. Finally, we determine the necessary physical 
quantities. We determine critical exponents below Tc in the present paper. 

Applying this method to the 0(N) invariant scalar model, the Euclidean Lagrangian 
density is given by 
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Here, the Ja are external source functions, and the index a runs from 1 to N. We assume that 
the coupling constant A is small, and therefore the perturbation theory at zero temperature 
is reliable. We first calculate the effective potential at an auxiliary large mass m = M ^ T 
at the one-loop level as 
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Here is a field expectation value. We leave only the finite-temperature part of the equation 
because we can ignore the zero temperature part due to the small coupling constant. We note 
that the daisy-resummation is not necessary because of the large mass. We then construct 
a non-perturbative evolution equation which connects the effective potential at an auxiliary 
large mass, ~ T^, and that of the true mass, = — /i^. Since we have constructed this 
for the 0(N) invariant scalar model in a previous work, & we present only the result: 
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This partial differential equation is solved with the initial conditions (2) numerically. 

We display the effective potential for N=4 around Tc in Fig. 1, and we find that the phase 
transition is of second order. The same behaviour is found for other values of N. This is 
consistent with other analyses using lattice field theory and renormalization group theory. & 
We find that the auxiliary-mass method satisfactorily deals with the problem of the infrared 
divergence. 
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Fig. 1. The effective potential obtained by the auxiliary-mass metliod (N = 4, A = 0.01). A 
second-order phase transition occurs at the critical temperature. Similar behaviour is observed 
at other values of N and A. 
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Fig. 2. Stable field expectation value (pc as a function of the temperature T (A = 0.01). 0c decreases 
monotonically and vanishes smoothly as the temperature increases. 




We next determine the critical exponents and study how well the auxiliary-mass method 
works. Since we have investigated this model above T^. previously, obtainig the critical 
exponents 7 and 5^ we investigate below and determine the critical exponent (3 here. 
The critical exponent, (3 relates an order parameter, (^c, to a reduced temperature, r = 



4 



as 
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The order parameter as a function of reduced temperature r is presented in Fig. 2 for 
N=4. Similar behaviour for other values of N is found. Since the order parameter (p^ 
vanishes smoothly at T^, we find that the phase transition is of second order. We next plot 
(pc as a function of r in Fig.3 for various N. These data appear linear with different gradients, 
corresponding to (] for each N. The exponent, is larger for larger values of N. 

We summarise the results of a present paper and a previous paper §^ in Table. I. The 
values of /?, 7 and 6 are much better than the Landau approximation and the dependence on 
N is close to the most reliable value (MRV). There are, however, slight differences between 
our results and the MRV, which are caused by an approximation in deriving Eq.3.0 

In conclusion, we have investigated the 0(N) invariant scalar model using the auxiliary- 
mass method and have obtained good results both qualitatively and quantitatively. These 
results suggest that the auxiliary-mass method is an effective tool at finite temperature. We 
were able to investigate not only second order phase transitions but also first order phase 
transitions since the finite-temperature field theory is based only on a statistical principle. 
We therefore believe that this is one of the most powerful methods to investigate a weak 
first order phase transition and models which have end-points : cubic anisotropy, the abelian 
Higgs model and the standard model.tfl 

Table I. The critical exponents (3, 7 and 6 obtained in a present paper and the previous paper. 
Those of Landau approximation (LA) and the most reliable values (MRV) are also summarised. 
We used the results of the six-loop approximation using the Pade-Borel resummation for the 



MRV here. 



/3 (LA,MRV) 



7 (LA,MRV) 6 (LA,MRV) 



N: 



N: 



N: 



N: 




0.39 (0.5, 0.327) 
0.41 (0.5, 0.348) 
0.44 (0.5, 0.366) 
0.45 (0.5, 0.382) 



1.37 (1, 1.239) 4.0 (3, 4.8) 

1.47 (1, 1.315) 4.2 (3, 4.8) 

1.60 (1, 1.386) 4.4 (3, 4.8) 

1.66 (1, 1.449) 4.4 (3, 4.8) 
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